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MODE  COMPETITION  IN  THE  QUASI-OPTICAL  GYROTRON 


I.  INTRODUCTION 

The  quasi-optical  gyrotron1-^  (QOG)  is  a  promising  millimeter  and 

q 

submillimeter  radiation  source  that  has  already  demonstrated  power 

levels  <  ISO  kV  vith  pulse  durations  up  to  13  ps  at  electronic 

efficiencies  <  14%.  Potential  applications,  ranging  from  heating  of 

fusion  plasmas  to  short  wavelength  radars,  require  stable  CV  operation  at 

high  power  levels  within  a  narrow  frequency  band.  Because  of  the  strong 

fields  developed  inside  the  cavity,  the  high  density  of  longitudinal 

resonator  modes,  and  the  relatively  long  time  of  operation,  nonlinear 

8  9 

interaction  among  many  modes  will  inevitably  occur.  '  In  fact,  it  is 
the  nonlinear  coupling  among  cavity  modes*1  that  will  determine  the 
existence  of  a  final  steady  state,  its  accessibility  from  the  initial 
small  signal  phase,  and  whether  it  is  a  single  or  a  multimode  state. 

In  short,  the  evolution  of  the  cavity  fields  can  be  described  as 
follows.  The  injected  electron  beam  initially  excites  all  the  cavity 
eigenmodes  characterized  by  a  start-up  current  below  the  beam  current. 
At  small  amplitude  each  mode  grows  exponentially  in  time  with  its  linear 
growth  rate,  unaffected  by  the  presence  of  other  modes.  Soon  one  or  more 
modes  reach  finite  amplitude  and  start  interacting  through  the  induced 
modifications  in  the  distribution  of  the  Kctron  beam.  Vhen  this 
happens  a  single  mode  may  nonlinearly  suppress  all  other  unstable  modes 
and  eventually  dominate.  Equally  well,  it  may  destabilize  modes  that  are 
linearly  stable.  In  some  cases  the  fastest  linearly  growing  mode  may  be 
overtaken  by  another,  linearly  slower  mode,  causing  mode  switching.  The 
final  steady  state,  if  one  exists,  may  involve  more  than  one  large 
amplitude  modes.  The  point  to  be  made  here  is  that  the  existence  of  a 
final  steady  state,  be  it  single  or  multimode,  is  determined  nonlinearly 
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and  cannot  be  predicted  from  the  linear  behavior. 

Finding  the  regimes  in  parameter  space  associated  with  the  desired 
type  of  final  equilibrium  proves  too  costly  for  direct  numerical 
simulation,  given  that  the  performance  of  the  QOG  depends  on  four 
independent  parameters.  It  should  be  stressed  that  the  final  equilibrium 
cannot  always  be  obtained  by  examining  the  linear  stability  of  every 
possible  single  mode  at  saturation. ^  Even  if  a  saturated  single  mode  is 
found  stable  to  small  perturbations  by  other  modes,  there  is  no  guarantee 
that  the  system  will  eventually  evolve  to  this  state.  In  general,  more 
than  one  stable  equilibrium  exists  for  a  given  set  of  modes,  however, 
only  one  of  them  is  accessible  from  the  appropriate  initial  conditions. 
In  the  case  that  the  final  steady  state  is  dominated  by  a  single  mode,  it 
is  not  necessarily  the  mode  with  the  largest  linear  growth  rate. 
Consequently,  all  participating  modes  must  be  treated  on  equal  basis 
until  (and  if)  one  dominates.  Mode  selection  and  accessibility  are 
inherently  nonlinear  processes. 

In  this  work  an  analytic  model  of  the  multimode  dynamics  is 
developed  by  expanding  the  nonlinear  current  in  powers  of  the  radiation 
amplitude.  A  set  of  equations  is  then  obtained  for  the  slow  time  scale 
evolution  of  the  wave  amplitudes  and  phases.  The  strength  of  the 
nonlinear  interactions  enters  through  coupling  coefficients  of  known 
analytic  dependence  on  the  gyrotron  parameters.  Numerical  simulations  of 
particle  trajectories  are  then  replaced  by  a  set  of  first  order  ordinary 
differential  equations,  one  for  each  participating  mode.  Not  only  is  the 
computation  time  reduced  by  orders  of  magnitude  but,  equally  important, 
the  final  equilibria  can  be  predicted  analytically  from  the  coupling 
coefficients.  In  this  paper  we  will  mainly  discuss  cases  with  up  to 
three  interacting  modes,  coupled  up  to  fifth  order  in  the  wave  amplitude. 
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The  coupled  equations  obtained  here  can  be  applied  to  the  description  of 
12 

phase-locked  operation  of  the  QOG  through  the  injection  of  a  small 

external  signal.  They  can  also  describe  phase  locking  through  the  use  of 

13 

a  prebunching  resonator,  via  certain  modifications  resulting  from  the 
use  of  a  prebunched  injected  electron  distribution  in  place  of  a  uniform 
one. 
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II.  GENERAL  FORMALISM 


The  configuration  for  the  QOG  is  shown  in  Fig.  1.  The  cavity 
fields  are  expressed  through  the  vector  potential  A(r,t), 


E  =  -  c  It  ’  B  =  7  x  A  ’  <D 

where  A  is  a  superposition  of  eigenmodes  of  the  Fabry-Perot  resonator, 

1  -  -iu  t 

A(r,t)  =  -  E  A  <t>  U  (r)  e  B  e  +  cc  ,  (2a) 

2  m  y 

Um(r)  =  Up  q  (r)  COskmX  *  ^ 

FmMm 

In  the  orthogonal  coordinates  (x,y,z)  with  x  along  the  resonator  axis  and 
z  along  the  electron  beam  axis,  the  transverse  profiles  upqCr)  are  given 
by  the  Hermite-Gaussian  functions 


P  Q 
Knf*m 


(r)  =  h  (y)  h  (z)  , 


m 


h„  <2>  -  "p  (iffy)  exp  (— 2~) 

*!n  Km  v  '  w  (x) 


(3) 


2  2  1/2 

w(x)  *  U  (1+x  /bm  )  is  the  radiation  spot  size,  V  is  the  radiation 
2 

waist,  b  «  k  V  /2  is  the  Rayleigh  length  and  H  (z)  are  the  Hermite 
mm  p 

polynomials.  The  wavenumber  km  is  equal  to  w„/c*  where  the  cavity 

frequencies  u  take  the  values 
m 


nc 

«  *  —  m- 

m  L  L 


2(p„+  q_  1 )  nc 

.  m  m 

-1  +  -  tan 


n 


(jt)]  =  -  <«  .  p..  ,„)  ,  (4) 


where  L  is  the  cavity  length  and  p_  and  q  are  integers;  the  notation  p. 

m  m  m 

q  is  used  for  labeling  frequencies,  u  *  o>  .  According  to  the  weak 
^m  m  mpq 
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coupling  approximation,  the  slowly  varying  amplitudes  Am(t)  are  complex 
quantities  with  magnitude  Ao(t)  and  phase  +m(t)> 

i*m(t) 

Vl>  =  •  •  <5) 

The  superposition  (2),  involving  various  frequencies  «m,  may  also  contain 
various  modal  structures  Pmqm  for  a  given  frequency  (i.e.  degeneracy). 


(a)  Particle  Equations 

Ve  use  the  guiding  center  description  for  the  transverse  particle 
coordinates 


x 


p  cos  6  , 


y  =  yg  +  p  Sin  9  , 


(6) 


PX  =  pgx  "  pi  sin  e  »  Py  *  pgy  +  pi  cos  9  ’ 

where  (x  ,y  )  and  (P^fP^)  denote  the  guiding  center  position  and 

momentum,  p^  is  the  magnitude  of  the  transverse  momentum,  8  is  the 

gyroangle,  p  =  P_l/0c  *s  the  Larmor  radius,  *  QQ/y  is  the  relativistic 

cyclotron  frequency,  where  SQ  «  |e|Bo/mc,  and  the  relativistic  factor  y  = 
2  2  1/2 

11  +  (p./mc)  +  (p  /me)  j  .  By  averaging  the  exact  Lorentz  force 

Z 

equations  in  the  vector  potential  representation  over  the  fast  (cyclotron) 
time  scale,  the  slow-time-scale  nonlinear  relativistic  equations  of  motion 
are  cast  in  the  form 


du. 

—  ■  E  “n  an  Ji(1V>  sin  <*n  *  V  cos  knxg  •  <7a> 

at  n 


de 

dt 


S  u 

o  r<  n 
—  +  L  —  a- 
r 


n  U1 


Jl<kn»> 

knp 


cos  (* 


+  ♦  )  cos  k  x 
n  n  g 


(7b) 
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In  Eqs.  (7)  time  has  been  normalized  to  w  a  (SQ/y  )  length  to 

kQ-  »  c/c»o,  an  *  lelAn/mc  is  the  normalized  radiation  amplitude,  u  is  the 

normalized  momentum  u  =  p/mc  *  yv/c  and  the  relativistic  factor  y  =  (1  + 
2  2  1/2 

u.  +  u„  )  .  The  angle  <|/  is  the  relative  phase  between  the  m-th  mode 

1  z  m 

and  the  particle,  =  0  -  »  t,  evolving  in  time  as 

mm 


dt 


-  6«  + 
n 


go 

En 

,T  3r 


U, 

n  1 


¥knp> 

knp 


cos  (* 


*  V  cos  knxg 


(8) 


where  Sa>  *  (a>  -  Q  /y)  «  1  is  the  zero-order  frequency  detuning.  The 
nnooo 

prime  (')  signifies  the  Bessel  function  derivative  in  respect  to  the 
argument.  The  evolution  of  y  is  found  combining  Eqs.  (7a)  to  (7b),  to 
obtain 


il  m 

dt 


Ji<k„»> sln  <*„  *  ♦„> cos  Vg  • 

n  v  ° 


(9) 


Two  approximations  are  implicit  in  Eqs.  (7)-(9).  The  guiding  center  drift 

has  been  ignored  and  the  axial  relativistic  momentum  uz  is  taken  as 

constant.  They  are  both  justified  for  an  «  1,  kV  »  1,  since  dx^/dt 

2 

dy  /dt  “  0(a  )  and  du  /dt  ~  0(a  /kV),  which  are  much  smaller  than  du/dt 

8  ^  *  n 

de/dt  *  0(an). 

(b)  Field  Equations 

Substituting  Eq.  (2)  into  Maxwell's  equations  and  ignoring  second 
2“  2  2~ 

order  derivatives  3  A  / 3 1  <<  «  A  ,  the  slow-time  evolution  of  the  vector 

n  n  n 

potentials  is  given  by 
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(10) 


The  second  term  in  the  right-hand  side  of  Eq.  (10)  is  zero  by  definition 

for  the  vacuum  eigenmodes  defined  in  Eq.  (2b),  thus,  all  the  information 

for  the  slov  evolution  of  the  fields  is  carried  in  the  nonlinear 

transverse  current  j^.  Equation  (10)  can  also  be  derived  from  the 

14 

general  expression  governing  coherent  wave  interaction, 


in  the  limit  of  standing  waves  inside  the  resonator  (zero  group  velocity, 

2  2 

Vg  =  0)  and  near  vacuum  dielectric  constant  e  (weak  beam  limit:  /«o  « 

1),  (de/dw)  a  2u>  . 

n  n 

Multiplying  both  sides  of  Eq.  (10)  by  exp(ia>nt)  and  taking  the  fast 

time  average  over  t  ~  2 n/to  isolates  the  nth  mode  on  the  left-hand  side. 

o  o 

Introducing  the  normalized  variables  used  in  Eqs.  (7),  and  letting  an  * 
|e|An/mc  =  an  exp  (i$n),  ve  have 


3a  i  4n|e | 

Vr>  if*  ~  1713  '  <12a) 

ot  /w  me 

o 

where 

1  _o  i«  t 

JNL(ri“n>  *  ~  J0  dt  JnL11’0  6  '  (12b) 

0 


Multiplying  both  sides  by  the  mode  profile  Un  and  integrating  over  the 
cavity  volume  lead  to 
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3i; 

i  2n}e |  f 

^  at 

2Q  J  n 

2  3 

<o  mcJV  J 

n 

o  n 

Vn  " 

J  d3r  Un2 

is  the 

phenomenological  term  «n/Qn  vas  added  to  account  for  resonator  losses  due 
to  boundary  effects,  where  the  cavity  Qn  is  defined  by 


d  ~  0 

—  a  *  =  -  it  a  2  .  (14) 

dt  n  Qn  n 

The  current  is  expressed  in  terms  of  the  nonlinear  distribution 
function  f^, 


-2n  .<*>  u. 

jNL(r,t)  =  “lel  nb<r>  Jq*10  J0duz  J0dUlUl  ~  COS0  fNL(r,U,t)  *  (15) 


Letting  fNL(r,u,t)  =  nQ  n^r^)  fNL(ultuz,z;t),  substituting  (15)  in  Eq. 
(13),  and  expressing  r^  m  (x,y)  and  »  (ux,uy)  in  the  guiding  center 
representation  Eq.  (6),  the  right-hand  side  of  Eq.  (15)  becomes 


2n  c  |e  |  n 


u  2  me3  V 
o  m 


7^  J^dxJ..d)rJ-.d2  "»(W  hPn‘V  V*’  cos 


V6  I0dUr  I0dUlUl  C0s6  Ji<kn0>  £NL<ui'u2'9’z;“n>  • 


where  again  fNL(n«n)  =  (1/*,)  J0T  dt  fNL^r,t^  exP(iwnt>-  Expression  (16) 
is  finally  recast  in  terms  of  the  injected  beam  current 


8 


h-  I' I  c(Jzono  J  dxg  I  %  WV  • 


(17) 


as 


/9  w i  .  - 

— - -  a  -  -  2ni  <  jMf(«  )  >  • 

V  2Q  J  n  JNL'  n' 

n 


(18) 


The  dimensionless,  volume  averaged  current  on  the  right-hand  side  of  Eq.  (18), 


Y  Tn  1 

<jNL(“n)>  *  - 

NL  n  u  c.V  T 
z  b  n  o 


rt 


0  i«„t 

dt  e 


dr  U(r)  jNL(r*t) 


lio_ 

"z  Vn  J 


dr  Un(r)  3»L<ri“n> 


(19a) 


is  given,  in  terms  of  fNL  from  Eq.  (16),  by 


I-  r 


Vn  uz  J° 


7  WV  hPn<J's>  “*  k"X 


n  g 


r®  p2n  .»  r®  ,  u, 

J_«dz  hqn(2)  J0de  J0dUz  VUlUl  Jl<knP)  “  cose  £nl<  We»,*V  • 


f-  3 

where  IQ  =  lell^/mc  is  the  dimensionless  beam  current  (Budker  parameter), 

and  a.  =  JJ  dx  dy  n.  (x  ,y  )  is  the  effective  beam  cross  section. 
d  g  g  d  g  g 


(19b) 
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III.  NONLINEAR  ELECTRON  DISTRIBUTION 


The  nonlinear  electron  distribution  in  phase  space  fM  (u. ,u  ,0,z,t) 
will  be  obtained  in  the  form  of  an  asymptotic  expansion  in  powers  of  the 
small  amplitudes  an  «  1.  Once  the  form  of  f^  in  the  multimode  field  (2) 
is  known,  the  nonlinear  transverse  current  j  ^  is  obtained  from  Eq.  (19). 
The  evolution  of  f„^  along  the  particle  trajectories  Eqs.  (7)  is  given  by 
the  collisionless  kinetic  equation 


3fNL 

3fNL 

de 

3fNL 

duz 

3fNL 

a. 

c 

1- 

!5 

sr*  “2 

4- 

3z 

dt 

39 

+  1 

dt 

3uz 

dt 

3u 

(20) 


By  formally  expanding  f  ^  and  ordering  in  powers  of  aJ ,  a  ~  0(aR)  «  1 


NL 


:(0) 


a  f 


(1) 


+  a 


j  f(j) 


one  obtains  a  hierarchy  of  equations 


L  f 


(j)  _ 


=  L,  f 


U-l) 


Df 


(j-1) 


where 


i  -  — 
o  "  3t 


+  u 


1.  +  s  L. 

3z  c  39 


(21) 


(22a) 


'1  =  Xa°  (21  Jl^knP> 


i*  3 

1  n 

cos  k  x  e  - 
n  g 


1  Jl(knp)  ,  4*n  8  1 

-  + - cos  ^  x  e  — 1 

3u .  2u .  k  p  n  g  3eJ 

1  in 

(22b) 


Operator  Lq  is  the  total  time  derivative  Lq  «  d/dt  along  the  unperturbed 
orbits  of  electrons  gyrating  in  the  external  magnetic  field  Bq 
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0'  -  0  -  ficT  ,  u±'  =  ,  Z1  =  2  —  *~T  t  *n'  -  *n  -  SwnT  ,(23) 

(the  characteristic  curves  of  LQf  =  0  in  the  theory  of  partial 
differential  equations)  where  0'  «  0(t'),  0  m  0(t)  and  t  =  t  -  t'. 
Replacing  the  exact  trajectories  in  the  right-hand  side  of  (21)  by  the 
unperturbed  orbits  (23),  (integration  along  the  characteristics),  it 
follows  that 

OD 

f(j)  =  J^dT  Df(j_1)  [u'(T),  0'(T),  Z'(T)  ]  ,  (24) 

We  exploit  the  fact  that  the  interaction  time  of  an  electron  inside  the 
cavity,  of  the  order  of  the  transit  time  t ^  =  W/0z,  is  much  shorter  than 
the  characteristic  growth  time  t  -1  =  d(lna  )/dt  for  the  fields.  Keeping 

cm  .-a 

aR  constant  during  the  integration  (24)  yields 


/•x  -  i(0  -  «  t.) 

f^  =  E  an  knX  e  J 

n  n  n  g 


(25) 


„n(2.Vi)  .“Vl  f  Jl(kn»)  8 


2i 


3u, 


JjO^p)  1  3  ' 

(knp)  2u1  30  , 


:(j-l) 


+  cc, 


where  Tj  *  *j  -  tj  ^  and  we  used 

-  Vj-1  -  ej  -  Vj  *  *Vj  • 

It  is  now  convenient  to  change  variables,  introducing  the  set  CQ,  Cj ,  ••• 
,  Cj,  defined  by 


u 


- - (t.+  x7  ♦  . . .  ♦  t.  ) 

W  YU  J 


(26) 
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where  dt.  =  -  (u  /yV)  dx.  and  C,  (x.*  0)  ■  C.  , .  Successive  iterations  of 
j  2  J  J  J  J-1 

Eq.  (25)  yield 


(27) 


where 


is  the  injected  beam  distribution,  and  is  written  as 


a 


i<e  -  vj* 


Vj> 


cos  k  x 
n  g 


W 


iv»< 


J1(knP)  3  Ji^knp)  1  3  ' 

2i  3^  (knp)  2u±  30  . 


(28) 


cc  . 


The  quantity  v  =  (Vy/u  )  5<o  ,  the  product  of  the  frequency  detuning  times 
the  transit  time  through  the  cavity  is  the  most  important  parameter  for 
the  mode  coupling  coefficients.  The  integrations,  Eqs.  (27)  and  (28),  are 
carried  out  in  Appendix  A.  The  expressions  for  f^\  j*l,2,3  are  given  by 
(A1)-(A3).  For  j  >  4  the  full  expressions  are  getting  rather  complicated. 
However,  one  may  exploit  the  existing  ordering  inside  f^  in  powers  of 
the  large  parameter  h  *  (kV/y)(u^/u2)  »  1.  It  is  then  easier  to  find  an 
approximation  for  f^  to  any  order  in  j,  Eq.  (A6),  by  keeping  the 
dominant  contribution.  Additional  simplifications  occur  during  the 
gyroangle  averaging  in  obtaining  the  nonlinear  current,  given  by  Eqs.  (A7) 
and  (A8). 
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IV.  MODE  COUPLING  EQUATIONS 


Substituting  expression  (A7)  for  the  nonlinear  current  inside  the 
right  hand  side  of  Eq.  (A8)  and  keeping  terms  up  to  fifth  order  in  complex 
amplitudes  am  =  am  exp(i$m),  the  general  mode  coupling  equations  assume 
the  form 


d^ 

dt 


2ni  y  I 

O  0 


uzo  °b  Vn 

(  Cm;n  am  +  Cklm;n  akalam  +  Cijklm;n  aiajakalam  )  * 


(29) 


A  mode  of  given  frequency  «n  interacts  through  all  possible  frequency 
combinations  satisfying  the  resonant  condition 


4"iJ....in  *  <“i  1  “j  1  •••  *  “„>  -  “n  *  0  ’  <30> 

A 

denoted  by  the  notation  {i,j,...,m)  implies  summation  over  all 

permutations  among  i,j,...,m  (except  n)  inside  Eq.  (29),  for  a  given  set 

of  resonant  frequencies  (30).  A  negative  index  -m  will  imply  -«m,  -$m, 

-\>  and  a  *  in  place  of  «  ,  ♦  ,  v  and  a„  respectively, 
m  n  m  m  m  n 

The  dominant  contribution  to  the  coupling  coefficients  is 

given  by 


2 

V*  S 


m 


r"  U12  (  y  N  r®  r2  df(0) 

i»  I0dUiT  ^  v"<z)  i-i*1  v"(Ci>  — 

zo  1 


Ola) 
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X 


,<3) 

■'klm 


;n  ■  (?)  1,4  W  “»  [‘(timin'  *  W)]  j’du.i-  (— )' 

^  zo 

T®  r2  ®  r^i  ®  r^2  ^  ^ 

J.d.Z  V"(2)  J>1  V«l>  —  Ld^  W  —  Jj«3  W  —  • 


(31b) 


5  2  5 

Cijklm;n  (?)  V  Sijklm;nexp [*  (^“ijklmjn*  +  *ijklmjn)]  J  du^-^—  (- — j 

^  20 

A2  vn<2>  JV  *i<ci>  “  J_S  yv  J-  £■••■  *zr  iS  w 


df 


(0) 


3u^  -®  3u^  -® 


du, 


where 


*  J1  <V>  hm«>  e 


-iv  C 
m 


(31c) 


(32) 


.  *  (♦•  ±  ±  • • •  ±  ♦  )  -  ♦  is  the  combined  slow  phase,  and  the 

lj . .  .m;n  i  j  ni  n 

nonlinear  filling  factor  is  given  by 


h.(yg)  hj(yg)...hn(yg)  cos  k^  cos  k^xg  ...  cos  kRxg.  (33) 

The  results  appearing  in  Eqs.  ( 29 )- ( 33 )  apply  for  any  spatial 
profile  and  velocity  distribution  of  the  electron  beam,  as  well  as  any 
radiation  modal  structure  in  the  open  resonator. 
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V.  THREE  MODE  INTERACTION 

Equations  (29)  are  quite  general,  involving  all  possible  resonant 
combinations  of  modes.  In  many  cases  it  has  been  observed  that  only  a 
small  number  of  modes  participates  in  the  final  equilibrium  with  finite 
amplitude.  In  this  section  the  three  mode  synchronous  interaction, 
satisfying  the  condition 

\  *  V1  *  2\  ’  or  “k  ♦  “l  *  2“m  •  (34> 

vill  be  considered  in  more  detail.  This  is  a  special  case  of  the  four 

mode  interaction,  with  «m  =  a>n  in  Eq.  (30).  Hence  the  frequencies  w^,  o>^ 
are  symmetric  around  «m,  <0^  =  «m  -  d«,  “  “m  +  dw,  where  d»  is  any 

frequency  interval  allowed  by  the  mode  spacing  in  the  cavity.  The 

possible  resonant  frequency  triads  entering  are 

Wm  '  “l  -  \ 

"  “m  -  “l  +  “1 

-  “m  -  \  +  \ 

-  to  -  to  +  to 

m  m  m 

and  all  3!  permutations  of  each  arrangement,  while  the  possible 
arrangements  for  are  given,  setting  i  *  k  and  j  -  1,  by 


9 


(35) 

9 


9 


'  “m  -  “1  +  “1  *  “k  +  \ 
-  wm  -  “k  +  “k  "  wm  +  wm 


“m  "  “l  '  “k  + 
“m  -  “l  -  “k  + 
wm  '  “l  '  “k  + 


“k 

w. 


“l 


w  -  to 


m 


m 


(36) 
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with  all  the  5!  permutations  of  each. 

There  has  been  a  tendency  recently  to  employ  the  following 

3 

renormalized  parameters  in  gyrotron  theory:  amplitude  Fn  =  an/Y03l0  = 
3  2 

E  /B  3  .  ,  detuning  A  =  (2/{5,  )$«  ,  interaction  length  u  = 

n  o  oi  n  10  n 

kV(0lo2/20zo)  and  current  I  =  IQ  (8n/Vo)(Q/r0Blo3)  where  VQ  *  (n/2)ko3LV2 

3  9 

and  I  *  Ib(Ams)  je  |/mec  x  3  10  .  Setting  =  w^,  and 

separating  real  and  imaginary  parts,  Eqs.  (29)  are  recast  in  the  form 


dF,  «, 


ur-UI-  J  n  fy  (  -  M\ 

-i  .  -i  fj  .  RTj  Fj  *  E  Pi  Fm2  «eGlm  .  F^  Fj  Rc («123  e'  1  * J 

at  m=  l 

E  E3  F,  Fj2  F.,2  ReDllm  .  E  F.2  F22  F3  R.  (Dml23  e'1*)  ,  (37.) 

1*1  m=l  m=i 


dF„  «, 


Wry  /  JA\ 

77  ‘  7T  F2  *  ReF2  F2  *  E  F2  V  ReG2„  *  F1  F2  F3  Re(G231  * ) 

dt  m=i 


3  3 


E  E,  F2  Fj2  F„,2  Rrt)2l.  *  E  Fm2  Fj  F2  F3  R.(d.231  eiF)  ,  (37b) 


1=1  m*l 


m=l 


dFj  w,  3 

dt  2Q 


ZT  F3  ■  ReF3  F3  *  E,  F3  Fn2  ReGlln  *  Fj2  Fj  Re(c312  **) 


m=l 

2  „  2 


*  E3  E3  F3  Fj2  F  2  R.D31m  .  E  Fm2  F  2  Fj  Re(om312  e’1*)  ,  (37c) 

1=1  m=l  m=l 


where  T  is  the  complex  linear  growth  rate.  The  evolution  of  the  slow  phase 
+  *  +2  +  +2  -  2$2  is  given  by 
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d* 

dt 


SS  +  Im 


(F1 


r3- 


2  r„ 


-  F12Im(G11+  G31-  2G21]  -  F32Im[G13+  G3;j-  2G23)  +  (G22“  Gj2-  G32] 


F3  F2 
F1 

3  3 


F,  F 


Io(G123  e  ”)  -  -h2-  Im(G312  e  “)  +  F1  F3  Im(G231  e'  ”) 


-  &  l,  Im(Dll«  +  D31m  -  2D21m)  *1 


^  “(w”) 


F1  F2 


-  i* 


Im(Dm312e  if)  *  F1  F3  Im(Dm231e 

(37d) 


where  5B  *  («o^  +  «3  -  Note  that  only  one  combined  slow  phase 
appears  for  each  triad. 

The  dependence  on  the  gyrotron  parameters  I,  y  and  A  is  contained 
in  the  coupling  coefficients.  It  will  be  shown  now  that  when  y  » 
yo0^o/2,  it  can  be  eliminated  as  an  independent  parameter  by  proper 
scaling.  The  computation  of  T,  G,  and  D,  carried  out  in  Appendix  B, 
yields 


rn  -  Vi  F<V  ’ 

-2 


nm 


nml 


J.S3  *  G<W  ’ 

IgS5  DCv^v.,^)  , 


(38) 


h»  mf  vl- 


.here  t  .  (ko«)2(Bi04/2Bz02)Jj'(k0r0)  end  Is  . 

I(0lo2/8Q)(koW)3(8lo/8zo)3IJ1'(kopo)]2.  The  nonlinear  filling  factors  Sr 
S3,  and  Sj  are  given  in  (Bll).  The  quantities  T,  G,  D  , given  for  a  cold 
beam  by  the  integrals  (B16)-(B18),  depend  on  the  desynchronism  parameter 
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\>m  alone,  Che. product  of  the  interaction  length  u  with  the  detuning  Affl, 


v  =  u  A 
m  m 


(39) 


This  suggests  the  following  scaling  transformation,  valid  when  u  =  «  and 

n  o 

Qn  =Q, 


I  *  2  IsQ 

P  *  l  F  ,  (40) 

T  =  t/2Q  , 

putting  the  mode  coupling  equations  in  the  final  reduced  form 

^7  *  F1  -  *  S3  Ej  F1F«2  B'G1«  *  S3F22f3  Re  (G123  •'  ") 

*  S5  Sj  FlFlV  ReDll«  *  S5  Re (Pnl23  «'**)}  '  <41» 


—  *  F2  *  I{SlBeF2F2  *  s3  m4  F2F«2  ReG2.  *  S3F1F2F3  Re(G123  e  **) 


dx 


*  S5  £,  iyiK*  R'°21*  *  S5  ^1%  Re(D„231  e“)}  •  <41b> 

^7  *  F3  -  I{SlReF3F3  *  S3  F3Fm2  R*G3„  *  Re(G312  •'  **) 


S5  £  £  VlV  ReR31m  *  S5  S  F.2F22Fi  »«(d„312  «-**) 
m*l  m= l 


(41c) 


It  follows,  from  (41),  that  the  normalized  start-up  current  is  1  = 

1/SlV 
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The  beam  parameters  I,  0ZQ,  the  interaction  length  u,  the 

cavity  Q  and  the  three  frequencies  combine  into  only  3  control  parameters: 
I,  Vj  and  \>2  (v3  *s  related  to  the  other  two  frequencies  through  the 

resonant  condition).  The  interaction  length  u  does  not  appear  explicity, 
but  only  implicitly  through  I  and  v.  For  simple  density  profiles,  such  as 
a  pencil  or  a  thin  annular  beam,  the  filling  factors  may  also  be  expressed 
as  powers  of  a  single  factor  s, 

Sj  *  s^+1  ,  s  *  cos  kQxg  h(yg)  ,  (42) 

with  x  and  y  the  effective  values  of  x  and  y  .  Thus,  they  could  also 
g  g  g  g 

be  absorbed  in  the  scaling  factor,  redefining 

V  -  s  l  , 

I'  =  s2  I  .  (43) 

Since  I  is  given  analytically  in  terms  of  the  gyrotron  parameters,  it 
suffices  to  compute  tables  of  T,  G  and  D  as  functions  of  once,  to 
completely  know  the  coupling  coefficients  for  any  combination  of  the 

parameters  I,  u,  0ZQ,  0^o  and  frequencies  wm.  Note  that  the  beam  spatial 
profile  enters  only  in  the  nonlinear  filling  factor  ;  in  case  of 

complicated  beam  profiles  Sj  can  be  computed  independently  from  Eq.  (33). 
The  discussion  so  far  has  been  limited  to  a  cold  beam.  In  case  of  thermal 

effects  the  general  expressions  (A7)  and  (A8)  are  applicable  and  the 

control  space  will  increase  by  two  parameters,  the  pitch  angle  spread  and 
the  energy  spread. 

Equations  (37)  and  their  reduced  form,  Eqs.  (41),  are  the  basic 
result  in  this  paper,  and  describe  three  mode  interaction  up  to  fifth 
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order  in  magnitude  in  the  QOG.  The  diagonal  terms  Gnn  describe  the 

effects  of  nonlinear  saturation  (or  self  excitation)  for  each  mode.  The 

cross-coupling  terms  Gnm  describe  a  phase- independe.  .  nonlinear  damping 

(positive  or  negative)  among  different  modes.  The  interaction  through  the 

G  ,  terms  does  involve  the  slow  phase  and  describes  mode  locking  effects, 
nml 

The  fifth  order  terms  D  ^  are  necessary  to  stabilize  the  system  in 
regimes  where  most  of  the  third  order  interactions  are  positive  and 
mutually  destabilizing.  They  are  also  required  to  account  for  potential 
amplitude  bistability,1^  i.e.,  the  existence  of  two  equilibria  of  the  same 
frequency  but  of  different  amplitude. 

Similar  results,  to  third  order  in  amplitude,  have  been  obtained 
for  the  conventional  gyrotron11,  using  a  somewhat  different  approach.  The 
present  derivation  of  the  coefficients  is  considerably  simpler,  since  it 
requires  a  (j+l)-fold  integration  for  the  j-th  order,  compared  to  the 
(2j+l)-fold  integration  required  in  the  single  particle  approach.11  Mode 
coupling  equations,  identical  in  structure,  first  appeared  in  the 
treatment  of  LASER  cavities1**. 
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VI.  RESULTS  AND  CONCLUSIONS 


The  normalized  linear  growth  rate  ReT(vn)  is  given  in  Fig.  2. 
Tables  of  the  third  order  coupling  coefficients  G(vn,\>m)  for  Gaussian 
profile  modes  and  a  cold  pencil  electron  beam  have  been  obtained  by 
numerical  evaluation  of  the  expressions  (B5)-(B7).  The  contour  plots  for 
the  real  and  imaginary  parts  of  Gnm,  G^  and  w^thin  the  regime  0  < 

vn,  vm  <  3  are  shown  in  Figs.  3a,  3b,  and  3c  respectively  (the 

coefficients  Gj12  are  Kiven  *n  terms  of  G^i  G3i2^vl,v2^  *  G123^v3’v2^ 

Dotted  lines  represent  regions  of  negative  values  for  G,  corresponding  to 

mutually  stabilizing  influence  among  interacting  modes  in  that  region. 

The  mode  coupling  is  destabilizing  in  the  regions  with  solid  lines, 

corresponding  to  positive  G.  Note  the  absence  of  symmetry  in  the  coupling 

coefficients,  where,  in  general,  G  t  G  .  Contour  plots  for  the  fifth 

nm  mn 

order  coupling  coefficients  D(vn*vm*vi>  as  functions  of  \>m,  for 
selected  values  of  \>n  appear  in  Figs.  4(a)-4(d). 

According  to  the  final  number  of  participating  equilibrium  modes, 
Eqs.  (41)  demonstrate  three  types  of  equilibria:  (a)  single  mode,  where 
one  mode  dominates  the  other  two  (b)  two  mode  equilibria,  where  one  mode 
is  of  negligible  amplitude  and  (c)  equilibria  among  three  modes  of 
comparable  amplitudes.  The  complete  set  of  equilibria  among  a  given  set 
of  frequencies  is  given  by  the  zeros  of  the  right-hand  side  of  (42). 
Since  Fn  are  defined  positive,  only  solutions  with  F^,  F2,  F^  >  0  are 

acceptable.  It  is  well  known  that  every  stable  equilibrium  of  a  nonlinear 
dissipative  system,  such  as  the  system  represented  by  Eqs.  (41),  is 
associated  with  a  region  of  initial  conditions  in  amplitude  space  (basin 
of  attraction)  that  eventually  fall  into  this  equilibrium.  Which 
equilibrium  the  system  will  choose  to  settle  in  cannot  be  analytically 
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predicted,  in  the  absence  of  invariants  of  the  motion.  This  can  only  be 
done  by  integration  of  the  set  of  ordinary  differential  equations  (41). 

The  simplest  case,  a  single  mode  equilibrium  for  the  n-th  mode,  is 
found  from 


?n  ( 


ReT  +  ReG 


nn 


ReD 


nnn  n 


0  , 


(44a) 


ImG 


nn 


ImD. 


nnn 


(44b) 


Equations  (44)  are  obtained  by  dropping  the  cross-coupling  terms  from 
(37);  the  combined  phase  #  is  irrelevant  in  that  case.  Acceptable,  non¬ 
trivial  solutions  Fn  >  0  for  the  steady-state  amplitude  should  always 

exist,  since  a  single  mode  saturates  through  particle  trapping.  The  fifth 

order  saturation  term  D  is  therefore  important  on  two  counts.  First, 

nnn 

it  is  necessary  for  saturation  in  regimes  where  the  third  order  GRn  is 

positive  and  causes  self-excitation  instead  of  suppression.  Second,  it  is 

required  to  account  for  amplitude  bistability,  when  two  acceptable 

equilibrium  values  Fn^  and  F^  exist  for  a  given  mode,  since  the  third 

1/2 

order  nonzero  solutions,  Fr  =  ±  (rn/Gnn)  ,  can  provide  at  most  one 

acceptable  solution  F  >  0.  At  steady-state  the  slow  phase  $n  varies  at  a 
constant  rate,  Eq.  (44b);  d$R/dt  expresses  the  nonlinear  frequency  shift 
from  the  linear  frequency  <i>n. 

Equations  (37)  are  used  to  provide  some  examples  of  mode 
competition  in  a  QOG  driven  by  a  cold  pencil  beam  in  a  typical  parameter 
regime,  with  beam  current  1^  =  13  A,  rQ  =  1.146,  3^O/0ZO  =  1*  V/A  =  5, 

cavity  Q  «  40,000  and  cavity  length  L  =  48cm,  corresponding  to  a  frequency 
separation  d«/«o  *  0.003.  In  Figs.  4(a)-4(c)  we  examine  mode  interaction 
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among  the  first  three  frequencies  <  “2  <  w3  a^ove  w0»  counting  only  the 
odd  resonator  modes  «n  separated  by  dw/w^.  =  0.006.  In  Fig.  5(a)  initially 
?2Q  »  F^q,  F^.  The  mode  F2  quickly  reaches  a  quasi-steady  state,  which, 
however,  is  unstable.  The  mode  F^,  which  initially  grows  much  slowly  than 
its  uncoupled  linear  rate  due  to  nonlinear  suppression  from  F^,  eventually 
overtakes  and  suppresses  F2  (mode  switching  in  the  cavity).  In  contrast, 
the  steady  state  reached  from  F^q  »  F2Q,  F.j0,  shown  in  Fig.  5(b),  is 
stable  since  a  large  amplitude  F^  suppresses  F2  and  F^  to  negligible 
amplitude.  However,  the  single  mode  equilibrium  of  F^  is  not  accessible 
from  the  usual  start-up  initial  conditions,  where  all  three  modes  have 
small  comparable  amplitudes  Flo  =  F2o  “  F3o  «  1.  Figure  5(c)  shows  that 
the  final  state  in  this  case  has  F^  as  the  dominant  mode.  Note  that  F2 
has  been  completely  suppressed  in  Fig.  5(c),  although  its  linear  growth 
rate  differs  from  F2  by  a  few  percent.  Since  F^  wins  the  competition  with 
Fj  and  F2,  we  next  examine  the  coupling  of  F^  with  the  next  two 
frequencies  F^  and  F^.  It  is  seen  from  Fig.  6  that  the  final  steady  state 
involves  three  modes  with  comparable  amplitudes.  Note  that  F^  has  been 
excited  despite  being  linearly  stable.  Figure  7  shows  the  mode 
competition  among  the  first  three  modes  in  a  cavity  of  reduced  length,  L  = 
24cm,  and  increased  frequency  separation  dco/ =  0.012  among  odd  modes. 

In  conclusion,  the  general  analytic  formalism  for  multimode 
interaction  in  the  QOG  has  been  developed  in  this  paper.  Applications  of 
the  coupling  equations  to  interactions  involving  three  main  modes  have 
demonstrated  the  ability  to  model  nonlinear  effects,  such  as  nonlinear 
destabilization,  suppression,  and  mode  switching.  A  new  scaling  has  been 
introduced,  during  the  derivation  of  the  coupling  equations,  that  reduces 
the  number  of  the  control  parameters.  The  ratio  of  operation  to  start-up 
current  I,  and  the  desynchronism  parameters  for  the  participating  modes 
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become  the  only  free  parameters.  The  inclusion  of  coupling  terms  up  to 
fifth  order  in  amplitude  was  shown  to  stabilize  the  system  in  the  hard 
excitation  regime  (where  the  third  order  coupling  coefficients  are 
positive.)  The  coupling  coefficients  up  to  fifth  order  have  been 
tabulated  for  a  cold  pencil  beam.  These  results  will  be  generalized  in  a 
future  paper  to  include  thermal  spreads  in  velocity,  and  arbitrary  beam 
current  profiles.  The  set  of  equations  will  also  be  expanded  to  include 
interactions  among  a  large  number  of  modes,  by  forming  all  the  possible 
resonant  frequency  combinations. 
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Appendix  A:  Computation  of  the  Nonlinear  Current 

Substituting  operator  (28)  inside  (27)  and  performing  the  successive 
iterations  yields 
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where  the  explicit  8,  t  dependence  is  contained  in  =  6  -  «nt  , 
r  ,  Jc  u ,  v  iv  C 

VC’U.L> - Jjf— J  cos  knx  hn(«  e  ,  x_n(^’ui)  •  Vc-“l> 

u  o 

zo 

Vc-“i>  -  v^i>  J1  ftr)  /  Ji fr1)  fir1)  -  ■  V<-V*- 


In  obtaining  (A1)-(A3)  we  used  the  invariance:  +  vfiC  =  const.,  along 

the  unperturbed  trajectories.  The  notation  (  -n,  m,  -1  )  inside  (A2)-(A3) 
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means  an  integral  similar  to  the  preceding  fully  written  term  except  that 
a  negative  index  -n  in  place  of  n  implies  -^n,  -\>n,  ^  aR*  in  place  of  \pn , 

V  *-„>  V  etc' 

The  dominant  contribution  in  each  integral  comes  from  the 

derivative  3/3u.  acting  on  v  , 
x  in 
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Therefore, 

the 

most  important 

terms 

inside  f^  are  those 

with  the 

(3/3u^)^  derivative.  Neglecting,  for  the  same  reason,  terms  like  3y/3u^  - 

1, 

3u.  /3u,  -  1  «  3v  /3u,  one  obtains  from  (A3) 
xx  m  x 


where  v  (£) 
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'if?) 
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s  =  cos  kx  h  (y  ) 
m  m  g  m  g 


One  can  show  by  induction,  proceeding  along  the  same  lines,  that 
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where  (...)  implies  all  the  permutations  among  the  m  ^ '  s. 

When  the  j-th  nonlinear  piece  of  is  substituted  inside  (15) 

i8  _  i  Q 

and  combined  with  cosG  ■  (l/2)(e  +  e  ),  the  resulting  expression 
contains  even  9-harmonics,  1,  e±^i6,  ...»  e±(j+*)i®f  for  j  0ddf  and  odd  0- 
harmonics,  e±i0,  e±3i0,  ...»  e±^+1^i®,  for  j  even.  Consequently,  coupling 

interactions  involving  even  number  of  modes  vanish  completely  during  the  0- 
integration.  The  first  nonlinear  correction  is  of  third  order  in  amplitude, 


the  second  of  fifth,  and  so  on.  It  should  be  emphasized  that  this  is  a 

consequence  of  assuming  a  uniform  in  0  initial  distribution  f^(u); 

however,  for  a  prebunched  distribution  of  the  form  f^(u,0)  =  fQ(u)  +  fj(u) 
16 

e  +  ...  +  cc,  both  even  and  odd  terms  survive  and  the  nonlinear  effects 


enter  to  second  order  in  amplitude.  Splitting  the  slowly  varying  fNL(«n) 
component  off  the  nonlinear  distribution,  by  multiplying  each  f^  with 
exp(i«nt)  and  averaging  over  the  fast  time,  substituting  ffjl/wn^  ^ns^e 
(19),  and  finally  averaging  over  the  gyroangle  0  and  over  the  resonator 
volume,  one  obtains 
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The  quantities  Aw  m  v„  *  ♦  and  S  are 

^2  *  *  *  ^  ®j  m^  m2  •  •  •  m^  { n  •  •  •  m^  }  n 

defined  in  a  similar  manner  as  in  equations  (30),  (32)  and  (33) 

respectively.  Of  all  the  possible  comoinations  [mi+n^tm^t. .  .±nu  J  only 
those  with  (j-l)/2  positive  and  (j+l)/2  negative  signs  have  survived  the 
©-integration  and  contribute  to  C^;  for  example,  the  term  t-l,-m,-n]  in 
f^3\  Eq.  (A5),  has  vanished  from  C^3^.  Also  note  that  the  complex 
conjugate  terms  in  (A1)-(A6)  drop  out  during  the  fast  time  averaging;  they 
contribute  to  the  jj4L^-wn^  component. 
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Appendix  B:  Computation  of  the  Three-Mode  Coupling  Coefficients 

The  three-mode  coupling  coefficients  in  Eqs.  (37)  are  computed  from 
the  general  expressions  in  (31),  keeping  the  frequency  combinations  (35) 
and  (36),  and  yielding, 
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A  cold  pencil  beam  is  employed  for  simplicity,  vhere  the  velocity 
distribution  f^  and  the  density  profile  n(x  ,y  )  are  given  by 

O  o 


f(0).  - -  S(ur  uiQ)  6(uz  -  uzQ)  ,  n<YV  “  5<Xg"?)  &(yg)  ‘ 

2nuLo 

The  calculations  are  simplified  by  the  following  observation  involving  the 
derivatives  of  Mu^-u^).  If  A(u),  B(u)  are  any  functions  of  u  satisfying 
dnA/dun  «  dnB/dun  for  every  h,  then 
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and  by  induction 
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Therefore,  moving  everything  but  the  exponentials  exp(i\>nt)  in  front  of  the 
velocity  integrals,  one  obtains 
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where 
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Ve  limit  the  mode  coupling  among  fundamental  Gaussian  modes  where  the 
transverse  profiles  Uqq  are 


uoo(yg’2)  =  w  ho(z)  » 


(B9) 


for  all  frequencies  a>n.  The  normalized  volume  in  that  case  is  given  by 
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For  the  pencil  beam  at  x  =  L/2,  and  according  to  Eq.  (4),  cos  k  x  =0  for 

g  n  g 

n  even  and  cos  k  x  >1  for  n  odd.  In  other  words,  the  electron  beam  passes 
n  g 

through  a  null  of  the  electric  field  for  modes  with  an  even  number  of  half 
vavelengths  betveen  the  mirrors,  or  through  a  maximum,  for  an  odd  number  of 
half  vavelengths.  The  j-th  order  nonlinear  filling  factor  S j ,  Eq.  (33),  is 
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then  reduced  to 
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Since  the  frequency  separation  among  the  coupled  frequencies  «n  is  very 

small,  it  is  important  only  when  the  detuning  S«n,  \>n  appears  explicitly. 

Thus,  one  may  set  kR  =  kQ,  letting  -  Ji'(kQP)  *  J0*  vn  =  v  = 

(n/2)kQ^LW^  and  Qr  =  Q  for  all  modes.  To  further  simplify  (B6)-(B10),  one 

may  notice  that,  according  to  (A4)  and  for  «n  =  «Q, 
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is  independent  of  frequency.  Since  the  S-function  causes  every  derivative 
to  be  taken  at  u^  «  u^Q,  one  may  employ  vr,  vm,  as  independent 
variables  and  use 


+ 


3  3  ^ 

»vm  +  3vl  J 


«  M  , 


to  show  that 


(B12) 


Fn  *  no  r(V  ’ 

(B13) 

-  3  * 

G„m  ■  "o  G<W  ■ 

(B14) 

°nlm  ‘  C  D<V  W  • 

(B15) 

36 


where 


n 


(B16) 


G<vv 


_®  ,«  p«  pC, 

-  Y,  I  J  dVndv«  j  «  vn  J  dCl  v-n 

{-n,m,-m} 


M 


(B17) 


D<vn'vl,',«) 


CD  CD  1 

I  lii 


dv  dv,dv 
n  l  d 


{n,-l,l,-m,m} 


_<DW  _CDW  _0D 


"  |C'dC4  vm  H  [  dC5  v-«  M  *(v  ''lo)  *(v  vno)  S(v  v.o)‘ 


(B18) 


It  is  clear  from  (B16)-(B18)  that  T,  G  and  D  depend  only  on  the  initial 
values  for  the  detuning  parameters  vnQ.  Combining  (B13)-(B15)  with  (B16)- 
(B18),  converting  the  amplitudes  an  into  the  normalized  Fn,  an  =  Fn 
(yo0^q  ),  and  absorbing  the  factors  (t0P^0  )  in  the  coupling  coefficients, 
finally  yields 
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and,  in  general,  for  the  j-th  order  coupling  coefficient, 
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puts  Eqs.  (B19)-(B21)  in  the  form  given  in  Eq.  (38), 
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Fig.  3  —  Contour  plots  for  the  real  parts  of  third  order  coupling  coefficients  (a)  (b)  Gm  and  (c)  G: 

Dashed  lines  signify  negative  values. 


Fig.  3  (Continued)  —  Contour  plots  for  the  real  parts  of  third  order  coupling  coefficients  (a)  G„  (b)  G,23  and 
(c)  G23l .  Dashed  lines  signify  negative  values. 
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2.5 


Fig.  4  (Continued)  - 
const,  plane  for  (a)  v 


1.0  1.5  2.0 


(d) 

-  Contour  plots  for  the  real  parts  of  the  fifth  order  coupling  coefficients  on  the  v, 
«  =  0.8  (b)  r„  =  1.9  (c)  vn  =  2.9.  The  contour  plots  in  (d)  are  drawn  for  =  •>, . 


0.2  0.4  0.6  0.8  1.0  1.2  1.4  1.6 


(C) 

Fig.  5  (Continued)  —  Mode  competition  in  the  QOG  for  the  first  three  modes  above  w„.  The  time  evolution  of 
the  amplitudes  F  in  (a)-(c)  corresponds  to  different  initial  conditions. 
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—  Mode  competition  among  the  third, 
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